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Abstract We adapt the method of Monte Carlo configuration interaction to 
calculate core-hole states and use this for the computation of X-ray emission 
and absorption values. We consider CO, CH 4 , NH 3 , H 2 O, HF, HCN, CH 3 OH, 
CH 3 F, HCl and NO using a 6-311G** basis. We also look at carbon monoxide 
with a stretched geometry and discuss the dependence of its results on the 
cutoff used. The Monte Carlo configuration interaction results are compared 
with EOM-CCSD values for X-ray emission and with experiment for X-ray 
absorption. Oscillator strengths are also computed and we quantify the mul¬ 
tireference nature of the wavefunctions to suggest when approaches based on 
a single reference would be expected to be successful. 

Keywords Monte Carlo configuration interaction • X-ray Emission • X-ray 
Absorption 


1 Introduction 

X-ray absorption can be used to experimentally study core-electron excita¬ 
tions, e.g., as has been applied to small organic molecules in Ref. [T]. While 
X-ray emission spectroscopy involves the initial ionization of a core electron 
followed by emission when the system adapts to remove the hole created in 
the core orbital. This experimental method has recently facilitated the investi¬ 
gation of dynamics in water [2] . More pertinent to this work it has previously 
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been used to probe the energy of relaxation of a valence electron back to the 
core in, e.g, simple alcohols and fluorine substituted methanes [1]. 

A method based on damped coupled cluster response has been created [S] 
to calculate X-ray absorption values and, for CCSD response, agrees well with 
experiment for neon, carbon monoxide and water. This approach would be 
expected to work well when the ground state is not considered multireference. 
Approaches based on density-functional theory (DFT) have also been created 
and shown to be accurate for small molecules, see, e.g., [HimiH] . These DFT 
methods can be applied to larger systems, however the functional used will 
affect the accuracy and current functionals are considered to not cope well 
with multireference systems. 

A successful computational approach to calculate the X-ray emission of 
many small molecules has been developed [Hiiin] using equation-of-motion cou¬ 
pled cluster singles and doubles (EOM-CCSD) [TT]. However for EOM-CCSD 
to be accurate the initial state should be able to be described well by CCSD, 
i.e., it should have a clearly dominant configuration when treated exactly in 
a given basis and therefore not be considered multireference. In the method 
of Refs. pilTU] a HF reference with a core hole is found using the maximum 
overlap method m then this is used for an EOM-CCSD calculation where the 
negative excitation energies are the emission values. Such an approach allows 
multiple emission values to be accessed in a calculation however there may be 
problems with the convergence [13] of the EOM-CCSD calculation and the ap¬ 
proach becomes intractable beyond reasonably-sized molecules. Furthermore 
if the full configuration interaction (FCI) core-hole wavefunction is deemed 
multireference then EOM-CCSD would be expected to neglect the static cor¬ 
relation of the core-hole wavefunction and so may have difficulties with excita¬ 
tion energies. Work on methods miM using time-dependent density-functional 
theory (TDDFT) offers the possibility of handling larger systems, but with a 
dependence on the approximations used and current functionals tend to have 
problems describing static correlation. For the further development of TDDFT 
approaches, in particular, the production of emission and absorption results for 
molecules of varying multireference character would therefore be useful. These 
data could also be used in improving the parameters in spin-component scaled 
configuration interaction with single substitutions and perturbative doubles 
SCS-CIS(D) [15] which also offers the possibility of emission calculations for 
larger molecules. 

Here we consider a complementary approach that is also limited to molecules 
that are not too large but should be able to deal with multireference situations 
and is not affected by convergence issues for a single emission calculation. To 
do this we adapt the method of Monte Carlo configuration interaction (MCCI) 
[niin] to describe core-hole wavefunctions. MCCI stochastically builds up a 
wavefunction with the aim of capturing many of the important aspects of 
the FCI wavefunction by accounting for both static and dynamic correlation 
to some degree, but using only a very small fraction of configurations. The 
method has been successfully applied to single point energies m. dissocia¬ 
tion energies [BUD], electronic excitations mm. ground-state mm and 
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excited potential curves m, multipole moments [H] and higher-order dipole 
properties up to the second hyperpolarizability m- 

We calculate the X-ray emission energies at equilibrium geometries for CO, 
CH 4 , NH 3 , H 2 O, HF, HCN, CH 3 OH, CH 3 F, HCl and NO. The emission en¬ 
ergy for CO at a stretched geometry of i? = 4 oq is also considered and we 
also look at the X-ray absorptions for the same set of molecules. We compare 
the emission energies with EOM-CCSD results of Ref. [5] when possible. These 
EOM-CCSD results have very good agreement with the available experimental 
studies. The absorption energies are compared in relation to available exper¬ 
imental results in the literature. The oscillator strength and multireference 
character of the states of interest are also computed and discussed. We note 
that this MCCI approach offers the possibility of multireference computational 
results for emission and absorption. We are not attempting to offer an improve¬ 
ment over EOM-CCSD for all systems and acknowledge that EOM-CCSD will 
be more accurate for systems that do not have significant multireference char¬ 
acter. However we hope that the results in this work will encourage tests, and 
possibly improvements, of EOM coupled cluster and TDDFT emission calcula¬ 
tions on more challenging multireference systems such as stretched geometries, 
nitric oxide and the carbon dimer. 


2 Methods 

MCCI [TOlfTT] randomly augments the configuration space by making single 
and double substitutions in the current selection of configuration state func¬ 
tions (CSFs) so that symmetry is preserved. By using configuration state func¬ 
tions the MCCI wavefunction is guaranteed to be a spin eigenfunction. The 
Hamiltonian matrix is then constructed using these configurations and diago¬ 
nalized. Any newly added configurations with an absolute coefficient, subject 
to appropriate normalization m, less than the cutoff (Cmin) are discarded and 
every ten iterations all configurations falling into this category are removed. 
After sixty iterations, the process continues until convergence in the energy, 
as described in Ref. m, is observed to 0.001 Hartree. The usual starting 
point is the configuration formed from the occupied Hartree-Eock molecular 
orbitals. The molecular orbitals and their required integrals are calculated 
using COLUMBUS p7] . 

Core-hole states could be calculated in MCCI by considering very high 
energy eigenvalues however for a stable calculation this would be expected to 
require all lower eigenvalues and so would not be feasible. Hence we extend 
the method to ground-state calculation restricted to a single occupied core 
orbital. 

Eor X-ray emission results, we initially perform a standard MCCI calcula¬ 
tion on the cation of the required symmetry. We then use MCCI to calculate 
the energy of the cation when the orbital containing the core hole is restricted 
to be singly occupied in all configurations. This is achieved by starting with a 
reference where the core orbital of interest is singly occupied then only allow- 
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ing substitutions that preserve this. One subtlety is that MCCI employs CSFs 
and uses the genealogical scheme m to ensure all orbital lists correspond to 
linearly independent CSFs. This means that the frozen single occupied orbital 
may be alpha spin in some lists and beta spin in others. As only non-frozen 
orbitals are available for substitution into existing configurations then for a 
randomly chosen configuration we check which spin does not have the frozen 
single occupied orbital and then allow the possibility of all but the double 
occupied frozen orbitals to be replaced in this spin. 

To calculate an X-ray absorption energy we begin with the neutral molecule 
then repeat the calculation for a core-hole state of the required symmetry with 
the lowest energy orbital singly occupied in all configurations. 

Below we summarize the use of MCCI for core-hole states starting with a 
reference consisting of a single occupied core Hartree-Fock molecular orbital. 

1. Create new configurations by random single and double substitutions in 
the current set of configurations so that symmetry and the frozen orbitals 
are preserved. 

2. Create the Hamiltonian and overlap matrices then diagonalize. 

3. Any new configurations with absolute coefficient less than Cmin are re¬ 
moved. 

4. Every ten iterations all configurations are considered as candidates for 
deletion. 

5. The procedure is repeated until the energy has converged. 

To calculate oscillator strengths between the two states of interest the 
following equation is employed 


fab = -AE\Dab\'^. 


( 1 ) 


Here 


Dab = {'I'a\r\'l'b) ■ 


( 2 ) 


We approximately quantify the multireference nature of the MCCI wave- 
functions by using the approach introduced in Ref. [23]. There 




(3) 


is calculated with an approximate normalization for configuration state func¬ 
tions such that X] = 1- Here a value of zero signifies that the wavefunction 
is single reference and one is approached as the system becomes more mul¬ 
tireference. Previous work [24] saw that the MR of an MCCI wavefunction 
for the strongly multireference chromium dimer when using a cc-pVTZ basis 
ranged from around 0.8 to almost 1 as the bond length was varied. In Ref. [25] 
the value for HF in a aug-cc-pVDZ basis was found to be 0.30 for an MCCI 
wavefunction suggesting that this system is amenable to being modelled using 
methods based on a single reference. 
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3 Results 

3.1 Emission Energies 

We first model the X-ray emission energy following the ionization of an electron 
from the lowest lying core orbital. If symmetry is used then both states are 
completely symmetric unless otherwise noted. We compare MCCI values using 
Cmin= 5 X 10“^ with the results of Ref. i for molecules that contain first row 
atoms and with one example (HCl) of a molecule containing a second row 
atom. We use the experimental geometry of the neutral molecule throughout 
except for methanol and CH 3 F where we optimize the geometry when using 
MP2 with cc-pVTZ. The calculations for CH 4 , NH 3 , H 2 O and HE used one 
frozen orbital while the other calculations used two except for HCl where five 
were employed. 


Table 1 MCCI emission energies and oscillator strengths at Cjjjjjj= 5 X 10“^ with a 6 - 
31IG** basis when using the lowest-lying core hole in the ionized molecule compared with 
experimental and EOM-CCSD results as listed in Ref. [2. 


Molecule 

MCCI (eV) 

f 


Fxp (eV) 

FOM-CCSD (eV) 

CO 

529.7 

7.1 X 

10-3 

- 

526.6 

CH 4 Ti ^ B 2 

277.3 

3.0 X 

10-3 

276.3 

276.2 

NH 3 

395.8 

4.3 X 

10-3 

395.1 

395.0 

H 2 O 

525.8 

4.3 X 

10-3 

525.1 

525.4 

HF 

674.7 

4.5 X 

10-3 

- 

674.5 

HCN 

393.9 

3.3 X 

10-3 

- 

393.1 

CH 3 OH 

529.7 

4.7 X 

10-3 

523.9 

522.2 

CH 3 F 

681.3 

3.0 X 

10-3 

675.6 

675.5 

CO {R = 4 uq) 

528.0 

1.5 X 

10-3 

- 

- 

HCl 

2821.0 

4.6 X 

10-3 

- 

2811.6 

NO 

544.0 

9.1 X 

10 -® 

- 

- 

C 2 (R = 1.25 A) 

287.8 

0 


- 

- 


The results for emission energies are displayed in table [T] and, except for 
CH 3 OH, CH 3 F and CO, are close to, but slightly higher than, the available 
EOM-CCSD results of Ref. [3] which themselves are in excellent agreement 
with experiment where available. This suggests that the core-hole state may 
not be described quite as well in this MCCI procedure as the cation. We note 
that these values for oscillator strengths are of a similar order of magnitude 
to those calculated with EOM-CCSD and a u6-311G** basis in Ref. [S] while 
HCN and CH 3 F used configuration interaction singles. The oscillator strengths 
demonstrate that the transitions are not forbidden within the dipole approx¬ 
imation except perhaps for NO. However for nitric oxide with a final state of 
Bi symmetry we find that / = 2.0 x 10“^ and the emission energy is 535.2 
eV while for a core hole in the second lowest orbital the MCCI emission value 
is 403.61 eV (/ = 4.0 x 10“"*) which is in good agreement with the experi¬ 
mental result m of X-ray lines around 403 to 402 eV assigned to a core hole 
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in N Is. For the latter core-hole state we find that MR = 0.87. This is an 
important result as it demonstrates that this approach can give good agree¬ 
ment with experiment when the multireference nature is very high. We also 
find that for the carbon dimer with a bond length of 1.25 angstrom that the 
multireference nature for the core-hole state is high at MR = 0.76 and the 
emission is 287.8eV. For the non-forbidden Ag i? 2 u transition we calculate 
/ = 4.7 X 10“^ and an emission energy of 289.4 eV. 

In table [5] we display the percentage error with the EOM-CCSD results. 
We see that there is very close agreement with the EOM-CCSD results, with 
the largest difference for methanol at 1.4%. HCl was considered in Ref. [9] 
using u6-311G** so we cannot compare directly and an experimental result is 
not available to our knowledge, but we note that their value was 2811.6 eV 
and that our result compared with this has an error of 0.3%. When using the 
cc-pCVDZ basis we calculate the emission as 2821.0 eV while the EOM-CCSD 
result [9] was 2805.9. A Hartree-Eock calculation with the Douglas-Kroll-Hess 
Hamiltonian in MOLPRO [29] suggests that in this basis the energy of the 
lowest energy Hartree-Eock orbital is reduced by 8.1 eV. This allows us to 
estimate the MCCI value when corrected for relativistic effects as 2829.1 eV. 
Eor the cc-pCVTZ basis the MCCI result is 2819.6 eV and the approximate 
correction for relativistic effects gives 2827.7 eV. 

As MCCI uses a random process to choose configurations we check that 
the results are sufficiently robust at this cutoff by repeating the calculations 
for water a total of ten times. We find that the mean emission energy is in 
agreement to one decimal place with the single calculation of table |T] at 525.8 
eV with a standard error of 0.0005 eV. 


Table 2 Percentage differences when compared with EOM-CCSD [9] when using MCCI at 
Cmin= 5 X 10“^ with a 6-31IG** basis when considering the lowest-lying core hole in the 
ionized molecule. 


Molecule 

Percentage Difference 

CO 

0.6% 

CH 4 

0.4% 

NH 3 

0.2% 

H 20 

0.08% 

HF 

0.04% 

HCN 

0 . 2 % 

CH 3 OH 

1.4% 

CH 3 F 

0.9% 


In tablejSjwe display the multireference values for the MCCI wavefunctions. 
When neither the cation nor the core-hole state is deemed multireference we 
display the molecular orbital transition when considering the most signifi¬ 
cant conhguration (|c| ^ 0.9) in each state. Eor the systems that we compare 
with the EOM-CCSD values at the neutral equilibrium geometry the core¬ 
hole cation would not be considered multireference except for perhaps carbon 
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monoxide. This suggests that the use of EOM-CCSD is indeed appropriate 
for these systems and even for carbon monoxide we note that the percentage 
difference is only 0.6% (table [2|) although in Ref. [28] there were two experi¬ 
mental emission values for carbon monoxide assigned to sigma orbitals at 522.3 
eV and 530.2eV, the EOM-CCSD result at 525.6eV lies between these values 
while MCCI is close to the higher value, however this was noted as being a 
very weak line. We note that the core-hole state is deemed multireference for 
NO suggesting that methods based on a single reference could perform poorly 
in this case. 

For carbon monoxide at a stretched geometry, both considered MCCI states 
are strongly multireference. We note that this stretched geometry results in 
a 1.7 eV change in the emission energy (tabled]). We investigate the effect of 
varying Cmin on the emission value for the stretched molecule. The FCI space is 
around 10® Slater determinants when symmetry is included while for the lowest 
cutoff considered we required 73883 CSFs for the cation and 115035 when the 
core hole is used. For Cmin= 5 x 10“'* 8702 and 10949 CSFs respectively were 
required. In Fig. |T] we see although the emission energy is non-variational it 
lowers with cutoff for the points considered. The plot suggests that for this 
challenging multireference system the results are still a little away from full 
convergence with respect to cutoff but we would not expect the emission energy 
to drop below around 527.6 eV. The emission energy reduces by around 0.35 
eV on lowering Cmin from the 5 x 10“'* value used for calculations in this paper 
to 1 X 10“'* and then by 0.04 eV to 527.64 eV on further reduction of the cutoff 
to 8x10“®. In Fig.|2]we see that for a system with low multireference character, 
hydrogen fluoride, there is again a decrease with cutoff but here the results 
seem much closer to convergence: the emission energy only reduces by 0.08 eV 
on lowering Cmin from 5 x 10“'* to 1 x 10“* and then by 0.004 eV on reducing 
Cmin to 8 X 10“®. 


Table 3 MCCI multireference character at Cj^jjj= 5 X 10 ^ with a 6-311G** basis for the 
ionized molecule with and without a hole in the lowest lying core orbital. 


Molecule 

Cation with core-hole MR 

Cation MR 

CO 

0.65 

0.35 

CH 4 1^2 lo-l 

0.43 

0.34 

NH 3 5a —>■ la 

0.33 

0.18 

H 2 O 3ai —>■ lai 

0.29 

0.15 

HF 3ai —>■ lai 

0.23 

0.11 

HCN 

0.45 

0.64 

CH 3 OH 

0.36 

0.75 

CH 3 F 9a —>■ la 

0.43 

0.33 

CO (R = 4 ao) 

0.80 

0.87 

HCl 5ai —>■ lai 

0.35 

0.20 

NO 

0.88 

0.56 

C 2 {R = 1.25 A) 

0.76 

0.72 
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Fig. 1 MCCI results for the emission energy of carbon monoxide {R = A ao) against 
on a logarithmic scale when using the 6-311G** basis. 
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Fig. 2 MCCI results for the emission energy of hydrogen fluoride against on a loga¬ 

rithmic scale when using the 6-311G** basis. 



Table 0] displays the percentage errors of our MCCI calculations and the 
EOM-CCSD calculations of Ref. [5] with the available experimental values 
listed in Ref. [3]. We see that EOM-CCSD is closer to the experimental results 
in the considered cases, however this is expected as this selection of molecules 
are not considered to have substantial multireference character for the core¬ 
hole state (table [3]). 

We use three systems as an example of the computational cost when using 
twelve processors for emission with increasing multireference character of the 
core-hole state. For HE the cation requires around 1 minute and 1320 CSFs 
while the core-hole state used around 5 minutes and 2724 CSFs. The CO 
cation required 13 minutes and 5838 CSFs and the core-hole state uses 53 
minutes and 11316 CSFs. For CO with a stretched geometry 55 minutes and 
8702 CSFs are needed for the cation while the core-hole state needed 1 hour 
and 38 minutes and 10949 CSFs. We note that in all three considered cases 
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Table 4 Percentage differences when compared with experiment when using MCCI at 
Cmin= 5 X 10“^ with a 6-311G** basis and for the EOM-CCSD results of Ref. [§|. 


Molecule 

MCCI 

EOM-CCSD 

CH 4 

0.35% 

0.04% 

NH 3 

0.18% 

0.03% 

H 20 

0.13% 

0.06% 

CH 3 OH 

1 . 10 % 

0.32% 

CH 3 F 

0.84% 

0 . 01 % 


the core-hole state is more challenging to compute and the cost increases with 
the multireference character. 


3.2 Absorption Energies 

We now consider X-ray excitation energies of an electron from the lowest lying 
core orbital in the same range of molecules rather than the emission energy. 
The results are presented in table [5) Unless otherwise stated, when symmetry 
is used we consider transitions between states classed as totally symmetric. For 
the Ai —>■ Ai transition in CH 4 , we find 289.0 eV compared with 287.1 eV for 
the experimental result [T]. The result stands out as the f value of 6 x 10“^^ 
indicates that this transition is forbidden within the dipole approximation. 
Hence we also calculate a core-hole molecule of B 2 symmetry. This gives 290.4 
eV and / = 3 x 10“^ while using the first Ai excited core-hole state gives 
290.8 and / = 3 x 10“^. We note that the experimental results range from 288 
eV to 290 eV for transitions assigned as to a t 2 orbital [T]. 

For water, the experimental absorption |30) for the first transition assigned 
to an Ai state is 534.0 eV which is also close to the MCCI calculation. Ab¬ 
sorption spectra for methanol have been calculated in Ref. m with the first 
peak at around 534 eV and a stronger absorption at about 537 eV which the 
MCCI result is close to. For CO the largest photoionization yield is between 
534 and 535 eV in Ref. [32]. The MCCI value is somewhat higher but this 
for an excitation of the same symmetry not an excitation to a tt orbital. For 
the excitation to Bi symmetry we find much better agreement as the absorp¬ 
tion energy is 535.6 eV with / = 3.5 x 10“^. For Nitric oxide Ref. [31] finds 
experimentally that absorption requires between around 532 eV and 534 eV 
for excitation to or ^Z\. These states are of A 2 symmetry when using 
C 2 V therefore agreeing with our result. For the B 2 B 2 transition we found 
/ = 2.9 X 10“® with an absorption energy of 542.9 eV. 

The damped coupled cluster linear response results of Ref. [5] for water and 
carbon monoxide are also in agreement with experiment with the exception of 
those from coupled cluster singles which are too high. The CCSD-NR result 
for water is 535.68 eV while for the CO excitation to a tt orbital it is 535.85 
eV. For water and carbon monoxide the multireference character is not high 
for the molecule (table jH]) suggesting that this approach would be expected to 
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be effective. For these absorption results, methods based on DFT with large 
basis sets have found 533.89 eV for water [7] while for CO the energies were 
534.21 eV [7], 533.0 eV [ 8 ] and, depending on the functional, 535.1 eV to 536.1 

eV [5]. 


,= 5 X 10-4 

a 6-311G** basis for the lowest energy core hole in the neutral molecule compared 
experimental results 


with 

with 


Molecule 

MCCI (eV) 

f 


Exp (eV) 

CO 

544.5 

1.4 X 

10-3 

538.9 

CH 4 

289.0 

5.9 X 

10-11 

- 

NHg 

402.8 

6.0 X 

10-3 

- 

H 20 

535.7 

1.2 X 

10-2 

534.0 

HF 

688.7 

1.9 X 

10-2 

- 

HCN 

405.2 

3.6 X 

10-3 

401.8 

CH 3 OH 

538.8 

9.3 X 

10-3 

534.1 

CHgF 

693.2 

1.4 X 

10-3 

- 

CO {R = 4 ao) 

531.2 

3.0 X 

10-2 

- 

HCl 

2832.0 

3.9 X 

10-3 

2823.9 

NO {B 2 A 2 ) 

534.1 

3.2 X 

10-2 

~ 532 to ~ 534 


Table [ 6 ] displays how all of the core-hole states for the neutral molecules 
appear to have multireference character. This continues in the core-hole state 
of Bi symmetry for carbon monoxide {MR = 0.71). The core-hole state 
also exhibits multireference character for the B 2 symmetry methane result 
MR = 0.66. Therefore a single calculation approach based on a single ref¬ 
erence would be expected to encounter difficulties. However earlier work [35] 
with unrestricted HF wavefunctions also achieved accurate results for CH 4 and 
remarked that the neglect of correlation cancels out to a large extent in this 
case when using the difference in energy between the neutral molecule with 
and without a core hole. 


Table 6 MCCI multireference character at Cj^jjj= 5 X 10 ^ with a 6-311G** basis for the 
neutral molecule with and without a hole in the lowest lying core orbital. 


Molecule 

MR 

core-hole MR 

CO 

0.47 

0.81 

CH 4 

0.37 

0.77 

NHg 

0.31 

0.74 

H 2 O 

0.32 

0.73 

HF 

0.26 

0.71 

HCN 

0.51 

0.86 

CHgOH 

0.45 

0.76 

CHgF 

0.41 

0.76 

CO (i? = 4 ao) 

0.89 

0.89 

HCl 

0.33 

0.78 

NO 

0.52 

0.75 
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As examples of the computational cost for absorption when using twelve 
processors we consider three systems of increasing multireference character. 
For HF the calculation for the molecule needed around 1 minute and used 
1339 CSFs while the core-hole state required 13 minutes and used 4164 CSFs. 
CO needed 5 minutes and 4576 CSFs for the molecule. The core-hole state 
needed 1 hour and 14 minutes and 10430 CSFs. For CO with a stretched 
geometry the molecule needed 45 minutes and 8836 CSFs while the core-hole 
state required 1 hour 47 minutes and 9639 CSFs. Similarly to the emission 
calculations the computational cost increased with multireference character 
and the core-hole states were more challenging. 


4 Summary 

We put forward a complementary approach to calculate X-ray emission and 
absorption energies for reasonably-sized molecules using Monte Carlo config¬ 
uration interaction (MCCI). This method should be able to cope sufficiently 
well whether the system is deemed to be well described by methods based on 
a single-reference or if multireference approaches are required. 

We saw that at equilibrium geometries the X-ray emission energies had very 
small percentage differences with the available EOM-CCSD results of Ref. |9]. 
When we quantified the multireference nature of the MCCI wavefunction we 
observed that the results suggested that the core-hole wavefunction tended not 
to be multireference in character and so EOM-CCSD would be expected to 
work very well for most of the systems. Nitric oxide was one of the exceptions 
to this where its core-hole state was deemed to be strongly multireference and 
the MCCI result for emission following the creation of a hole in the second 
lowest energy orbital compared well with experiment. This suggests that sim¬ 
ilar open shell systems may pose difficulties for emission calculations when 
using approaches built around a single reference. We also considered carbon 
monoxide at a stretched geometry and saw that the system would be consid¬ 
ered multireference with an accompanying change in the X-ray emission of 1.7 
eV. 

We also looked at the X-ray absorption of the molecules and compared the 
MCCI results with experimental data when available. For methane we found 
reasonably good agreement with experiment for the excitation of an electron 
from the lowest lying core orbital. The results with water, ethanol, hydro¬ 
gen cyanide and nitric oxide also fitted in with known experimental values. 
The value for hydrogen chloride was about 8 eV higher than experiment. The 
largest absorption energy in carbon monoxide was higher than experiment 
but for excitation to Bi symmetry we found much better agreement with ex¬ 
periment. Interestingly the multireference character of the core-hole MCCI 
wavefunction was fairly large implying that methods based around the unre¬ 
laxed core-hole single-reference may encounter difficulties for these absorption 
calculations. 
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This approach can be straightforwardly extended to consider holes in or¬ 
bitals that are not the lowest in energy and we have illustrated this on nitric 
oxide. When each wavefunction has only one significant configuration then we 
can label the transition using two molecular orbitals, however when dealing 
with multiconfigurational wavefunctions although we choose the core-hole or¬ 
bital, it is not trivial, or perhaps possible, to label the transition in terms of 
a single excitation using molecular orbitals. The use of natural transition or¬ 
bitals [36] , as used for wavepackets created by X-rays m, or natural transition 
geminals [38] may allow this transition to be assigned a compact description. 

These calculations of X-ray emission and absorption for reasonably sized 
molecules with strong multireference character should provide useful data 
for improving approximations in methods for larger systems such as time- 
dependent density-functional theory. 
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